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We study numerically the phase diagram and the response under a driving force of the phase field crystal 
model for pinned lattice systems introduced recently for both one and two dimensional systems. The model de- 
scribes the lattice system as a continuous density field in the presence of a periodic pinning potential, allowing 
for both elastic and plastic deformations of the lattice. We first present results for phase diagrams of the model 
in the absence of a driving force. The nonlinear response to a driving force on an initially pinned commensu- 
rate phase is then studied via overdamped dynamic equations of motion for different values of mismatch and 
pinning strengths. For large pinning strength the driven depinning transitions are continuous, and the sliding 
velocity varies with the force from the threshold with power-law exponents in agreement with analytical pre- 
dictions. Transverse depinning transitions in the moving state are also found in two dimensions. Surprisingly, 
for sufficiently weak pinning potential we find a discontinuous depinning transition with hysteresis even in one 
dimension under overdamped dynamics. We also characterize structural changes of the system in some detail 
close to the depinning transition. 

PACS numbers: 64.60.Cn Order-disorder transformations; statistical mechanics of model systems, 64.70.Rh Commensurate- 
incommensurate transitions. 68.43.De Statistical meclianics of adsorbates, 05. 40. -a Fluctuation phenomena, random pro- 
cesses, noise, and Brownian motion 



I. INTRODUCTION 

There exist many systems in nature with two or more com- 
peting length scales, which often leads to the appearance of 
spatially modulated structures. Such systems may exhibit 
both commensurate (C) and incommensurate (I) phases lU 0] 
characterized by differences in the spatial ordering of the sys- 
tem. Important examples include spin density waves iH 01, 
charge density waves vortex lattices in superconducting 
films with pinning centers and weakly adsorbed monolay- 
ers i?! 01 on a substrate. The emerging structures are char- 
acterized by an order parameter (e.g. charge, spin or particle 
density) that is modulated in space with a given wave vector 
q. In particular, for 2D adsorbate systems, there is compe- 
tition between the commensurate state which is favored by a 
strong periodic pinning potential and the cost of the elastic en- 
ergy depending on the mismatch between the intrinsic lattice 
constant a of the overlayer, and the period b of the pinning 
potential. 

While the static properties of C and I structures have been 
extensively characterized iH 01 much less is known about 
their dynamics. A particularly interesting case arises when an 
initially pinned phase is subjected to an external driving force 
/. The resulting nonlinear response is relevant for a variety 
of different physical systems which are accessible experimen- 
tally. A driven atomic monolayer in a periodic pinning po- 
tential is an interesting realization of such nonlinear behavior 
Jgtl, which is directly relevant for experiments on sliding fric- 
tion between two surfaces with a lubricant lUOll a nd between 
adsorbed layers and an oscillating substrate ill ll [l2ll . Other 
systems of great interest are driven charge density waves 
[Tsll in which commensurability and impurity pin- 



ning often compete lfT6ll . and superconductor vortex arrays in 
which different commensurability and pinning behaviors have 
been experimentally observed lfl7l[l8ill9ll . including periodic 
and asymmetric potentials. 

For a sufficiently large pinning potential, the phase may re- 
main pinned for small forces if there are no thermal fluctua- 
tions present. This means that at zero temperature there is a 
finite critical force fc above which the system starts moving. 
For many systems, it is found that just above the threshold 
fc, the drift velocity Vd shows a power-law dependence with 
respect to the force / 

«d«(/-/c)^. (1) 

When this behavior is regarded as a dynamical critical phe- 
nomenon, the power-law exponent ^ of the corresponding 
driven depinning transition can be argued to result from the 
scaling behavior of the system near the threshold ITsIl with 
corresponding divergent time and length scales and universal 
behavior In general, however, the observed value of ( may 
depend on the system and its dimensionality. For a pure elastic 
medium with quenched randomness there appears to be a uni- 
versal value which depends on the dimensionality of the sys- 
tem, provided inertial effects are negligible ll20ll2lll . For the 
case of an initially commensurate phase in a periodic pinning 
potential without disorder, a power-law exponent ( = 1/2 is 
expected independent of the dimension as the threshold be- 
havior can be understood from the point of view of single par- 
ticle behavior [Tsi I22I1 . In the limit of a large force 
(/ ^ fc) / fc ^ 1 the system is moving and the corresponding 
relationship between the driving force and the sliding velocity 
defines the sliding friction coefficient 

Vs = f/vd- (2) 
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The simple Frenkel-Kontorova (FK) model 11231 12411 ex- 
tended to two dimensions and other similar elastic models 
have been used to study driven depinning transitions and the 
sliding friction of adsorbed monolayers |25l]. Although 
these models take into account topological defects in the form 
of domain walls they leave out plastic deformations of the 
layer due to other defects such as dislocations. These defects 
are particularly important when the CI transition occurs be- 
tween two different crystal structures, or in presence of ther- 
mal fluctuations or quenched disorder, and should be taken 
into account for a more realistic description of the system. 
Such defects can be automatically included in a full micro- 
scopic model involving interacting atoms in the presence of 
a substrate potential using more realistic interaction poten- 
tials. However, the full complexities of the microscopic model 
severely limit the system sizes that can be studied numeri- 
cally, even when simple Lennard- Jones potentials are used to 
describe the interactions 1126112711. 



When the driving depinning transition is discontinuous, 
hysteresis effects can occur which result in two different crit- 
ical forces, > Z^'^, corresponding to the threshold val- 
ues for increasing the force from zero and decreasing the 
force from a large value, respectively. A fundamental issue 
in modeling such systems is the origin of the hysteresis. It 
is well known that hysteresis can occur in underdamped sys- 
tems, where inertia effects are present. However, molecular 
dynamics simulations of a 2D model of an adsorbed layer with 
Lennard- Jones interacting potential for increasing values of 
the damping coefficient (microscopic friction) and analytical 
arguments suggested that hysteresis should remain in the over- 
damped limit 13, IS- On the other hand, results for the pure 
elastic FK model shows that although the hysteresis behavior 
is similar to Lennard-Jones model for weak damping, it dis- 
appears in the overdamped limit ll25ll . The different behaviors 
could be due to the absence of some defects generated during 
the depinning transitions, which are allowed in the Lennard- 
Jones model but not in the elastic FK model. In fact, hys- 
teresis can be argued to arise from topological defects in the 
lattice such as dislocations in two-dimensional systems lfl4ll 
even in the absence of inertial effects. One thus expects that 
overdamped dynamics should be able to described the hys- 
teresis behavior in two dimensions provided the model incor- 
porates both elastic and plastic deformations. For charge den- 
sity waves, where the pinning potential is disordered, a field 
theoretical model has been introduced which allows for dis- 
locations as well as thermal fluctuations ll29ll through ampli- 
tude and phase fluctuations, and shows both elastic and hys- 
teretic behavior in agreement with experiments jsoll . In ab- 
sence of disorder, however, the possible hysteresis behavior in 
such models has not been investigated. 

Recently a phase field crystal (PFC) model was introduced 
lIsTl l32i [33I1 that allows for both elastic and plastic deforma- 
tions in the solid phase. In this formulation a free energy func- 
tional is introduced which depends on the field il){r,t) that 
corresponds to the particle number density averaged over mi- 
croscopic times scales. The free energy is minimized when 
■0 is spatially periodic {i.e., crystalline) in the solid phase 
and constant in the liquid phase. By incorporating phenom- 



ena on atomic length scales the model naturally includes elas- 
tic and plastic deformations, multiple crystal orientations and 
anisotropic structures in a manner similar to other microscopic 
approaches such as molecular dynamics. However, the PFC 
model describes the density on a diffusive and not the real mi- 
croscopic times scales. It is therefore computationally much 
more efficient. 

In our previous works ll34l Issll we demonstrated how the 
influence of an external periodic pinning potential can be in- 
corporated in the PFC model. Such a model provides a con- 
tinuum description of pinned lattice systems. The pinning po- 
tential is chosen such that it allows the occurrence of both C 
and I phases as ground states of different symmetries in the 
model. In Ref. Isill part of the phase diagram as function 
of pinning strength and lattice mismatch between the pinning 
potential and the PFC was mapped out. Numerical minimiza- 
tion was used to find the minimum free-energy configurations 
and provide details on the topological defects in the bound- 
ary region. In particular, we found that the transition from the 
I to the C phase remains discontinuous for all values of the 
mismatch studied in Ref. ll34ll . We also performed a detailed 
Voronoi analysis of the defects throughout the transition re- 
gion. In Ref. issll the equilibration method was improved and 
the range of mismatches extended to include both positive and 
negative mismatches. 

In the present work we focus on the case where the PFC un- 
der an external periodic potential without disorder introduced 
in Ref. ||34ll is driven by an external force in the absence of 
thermal fluctuations. To this end, we first present improved 
detailed phase diagrams of the model both in ID and 2D. The 
main focus of the present work is on the influence of an exter- 
nal driving force on the pinned C phase, which we study for 
different values of mismatch and pinning strengths for ID and 
2D systems. As expected, due to the competition between the 
pinning potential and the driving force there is a depinning 
transition at fc for a finite driving force /. We demonstrate 
that within a certain range of parameters the depinning tran- 
sitions are continuous, and find that both in ID and 2D the 
corresponding power-law exponent is = 0.5 in agreement 
with the expected value iO [l4l flSl 12211 . We also characterize 
structural changes of the system close to the depinning tran- 
sition. For large pinning strength transverse depinning tran- 
sitions in the moving state are also found. Surprisingly, for 
sufficiently weak pinning potential we find a discontinuous 
depinning transition with hysteresis even in one dimension al- 
though overdamped dynamical equations are used. 



II. THE PHASE FIELD CRYSTAL UNDER A PERIODIC 
POTENTIAL AND A DRIVING FORCE 

For the phase field crystal in the presence of pinning poten- 
tial Ii31, 32, 33^^J]^ the free energy functional can be written 
in dimensionless form as 



F 



dx 



V{x)il;{x) 



(3) 
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where r is a temperature dependent quantity and V{x) is an 
external potential which represents the effect of the substrate. 
This model can be derived directly from the classical density 
functional theory llssll of freezing by expanding around the 
properties of a liquid in coexistence with a solid phase. More 
specifically it can be shown ll36ll that r is proportional to the 
difference between the isothermal compressibility of the liq- 
uid and the elastic energy of crystalline phase. Furthermore 
the length scales in this model have been scaled by the nearest 
neighbour distance in the coexisting liquid state. 

In the absence of the pinning potential the equilibrium 
minimum energy configuration of the system depends on 
the parameter r and the average density "tp = ^ J dx'ip{x) 
lf32ll . where Vd is the system volume in d dimensions. In 
2D, the solid phase corresponds to a triangular lattice. The 
length scale chosen here corresponds to k^^, where fco = 
2TT/{atV3/2) = 1, with at as the lattice constant of the in- 
trinsic triangular lattice. When the external pinning potential 
is present, the competition between the length scales associ- 
ated with the intrinsic ordering and the pinning potential can 
lead to complicated phases depending on the parameters cho- 
sen in the energy functional P4i , i35[| . We choose an external 
pinning potential of a simple periodic form V = Vq cos{ksx) 
in ID, and V = Vo[cos(fcsx) + cos(fcsy)] in 2D. The wave 
vector kg is related to the periodicity of the pinning potential 
tts, such that ks ~ 27r/as. 

We define the relative mismatch between the external 
potential and the PFC as 



Sm = (1 - kg) 



(4) 



The response of the system to a driving force can be obtained 
by includin g a convective derivative, / • Wtp, to the original 
PFC model [331 ■ Thus, the dynamical equation of motion for 
the phase field is given by 

OT dip 

= ((r + (1 + V^)') V.' + i'^ + V^ 

+/ • Vt/;. (5) 

For the forces considered in this work (typically / = fx), 
the convective term does not change the average value of the 
density field. 

In contrast to the usual classical microscopic characteriza- 
tion of particle positions and velocities, the measurement of an 
average drift velocity in response to an external driving force 
/ requires some discussion. In the PFC model the maxima of 
the density field that define the lattice structure cannot always 
be interpreted as individual particles, since vacancies may be 
present in the system. The conservation law in the model con- 
cerns the local density field, not the number of maxima in the 
field. This becomes evident in the driven PFC model, where 
the motion of the density field close to depinning may be more 
akin to flow in a continuous medium than the motion of dis- 
crete particle-like objects. Thus, defining the drift velocity 
in terms of the density field maxima is computationally more 
demanding to implement. We have found that measuring the 



drift velocity Vd from the rate of change of the gradient of the 
density field gives consistent results, in absence of thermal 
fluctuations. We have used the following definition 



vd^{m/dt\),/{\d^/dx\),)^, 



(6) 



where the subscripts x and t in the brackets denote averag- 
ing over space and time, respectively. Alternative definitions 
were also considered, but this particular form proved the most 
statistically accurate. 

Although the definition of the drift velocity Vd according to 
Eq. (|6]l can be used to determine the velocity response along 
the direction of the driving force, it is not particularly useful 
in the study of the response in the transverse direction since 
it is not a vector quantity. In order to study the transverse re- 
sponse it is more convenient to determine the average velocity 
directly from the positions of the local peaks in ip{x). This 
requires locating such peaks as a function of time during the 
numerical simulation. We have developed a computational 
method which determines the location and velocity of each 
individual peak in presence of the external force and thermal 
fluctuations fyf\\. The method to locate the peaks is based on 
a particle location algorithm used in digital image processing 
OSll . The drift velocity for the lattice of density peaks is ob- 
tained from the peak velocities Vi as 



Vp 




(7) 



where Np is the number of peaks. We find that the definitions 
of the velocity from Eqs. ^ and Q give consistent results 
for the longitudinal depinning in absence of thermal fluctua- 
tions. However, in presence of thermal fluctuations only the 
definition from the peak velocities O is able to separate the 
contribution to the drift velocity due to the driving force from 
thermal noise contributions. 



III. RESULTS 

In this section, we present results obtained for the static and 
dynamic properties of the PFC model. Numerically, we study 
the system properties by integrating Eq. (|5]) using a simple 
Euler algorithm and the time derivatives are approximated by 
a forward finite difference with the time step dt = 0.005 (the 
time scale corresponds to the diffusion time over the length 
scale k^^). For the ID case, the density field was discretized 
on a uniform grid with dx = tt/A, while for the 2D case we 
used a square uniform grid with dx = dy = -k / A. The lapla- 
cians are evaluated in ID using a central difference, while 
for the 2D the 'spherical laplacian' is used ll32l l39ll . In both 
the equilibrium and driven situations, fully periodic bound- 
ary conditions have been used. Note that for the conserved 
time-dependent Ginzburg-Landau (TDGL) equation (|5]l, the 
addition of the driving force of the form f{dip/dx) preserves 
the local conservation of ip under periodic boundary condi- 
tions. For the study of the static properties / is set to zero. 
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For a given value of the mismatch, the pinning strength Vq 
is increased from zero to a maximum in steps of cIVq and 
then decreased back to zero. Each time the pinning strength 
is changed we allow the system to equilibrate. The final 
state corresponds to a configuration that minimized the energy 
functional. 

For the study of the influence of a driving force we choose 
the mismatch and pinning strength such that the system is ini- 
tially in a commensurate state. The driving force / is then 
increased from zero to a maximum value | / 1 > fc and then 
decreased back to zero. In the case where the depinning tran- 
sition at fc is continuous, we determine the corresponding de- 
pinning exponent ( in the limit |/| — > fc. Unless specified the 
force is applied in the x direction, i.e., / ~ fx. 

A. Phase Diagram and Nonlinear Response for ID System 

1. Equilibrium properties 

The static properties of the pinned PFC in ID (/ = 0) are 
presented as a phase diagram in the Vq — Sm plane shown in 
Fig. [T] For comparison, we include in Fig. [T]the phase bound- 
aries obtained analytically and numerically from from Eq. |5] 
The analytical phase boundary was obtained by minimizing 
the free energy F[ip{x)], expanding the density field as 

tp{x) = Ai cos{x) + A2 cos{ksx) + A3 cos((2 — ks)x), (8) 

where the last term accounts for the distortion of the lattice. 
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FIG. 1: Phase boundary between commensurate (C) and incommen- 
surate (I) phases for the ID pinned PFC model, calculated numer- 
ically (continuous line with triangles) and analytically (continuous 
line). 



Next we investigate the influence of an external force on the 
ID commensurate phase. For this purpose, the parameters are 
chosen such that r = — 1/4, and ij; = 0. Depending on the 
values of the mismatch and pinning strength different behav- 
ior is found when the driving force is added. Several values of 



FIG. 2: (a) Discontinuous depinning of the commensurate phase for 
relatively low pinning strength (5ni = 0.3125, Vo = 0.11). (b) De- 
pendence of A/c on Vo for 5m ~ 0.3125. 



mismatch between 0.3125 and —0.50 were investigated. For 
i5m > —0.3 two types of depinning are present, discontinu- 
ous and continuous. For values of the pinning strength close 
to the I-C phase transition, the depinning transition is discon- 
tinuous. The dependence of the drift velocity as a function 
of the driving force exhibits a hysteresis (see Fig. |2(a)| i. The 
gap A/c = /™ — fc'^ decreases when the pinning strength 
increases (Fig. |2(b)[ ), which indicates that the transition be- 
comes continuous for large enough Vq. 

Note that for S^^ < —0.3 the depinning transition of a 
C pinned phase is continuous for all values of the pinning 
strength, while for ^ —0.3 only for values of pinning 
strength above a certain threshold. In Fig. [3] we show the be- 
havior of a continuous depinning transition for 5^ = 0.3125, 
Vq = 0.11. The dependence of the drift velocity on the force 
follows a power law Vd oc {f — fc)'', as can be seen in Fig. 
[3] The exponent C does not depend of the pinning strength 
and it is equal to 0.50 ± 0.03 in all cases studied here. This 
result can be understood as follows. When the pinning poten- 
tial is periodic and large in magnitude the neighboring phases 
are weakly coupled and the system should behave as a single 
particle in a periodic potential ifisl [Til [Tsi I22I1 . This effec- 
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FIG. 3: Continuous depinning transition for relatively high pinning 
strength (Sm ~ 0.3125 and Vb = 0.250). The arrow marks the value 
of the critical depinning force. The triangles represent the numerical 
data while the continuous line is a fit to (/ — /c)^. The best fit is 
obtained for the exponent — 0.50 ± 0.03. 



c(2 X 2) phase (Fig. |4(d)| i in which every second site of the 
lattice of the pinning potential corresponds to the maximum in 
the phase field This state is favored for mismatch values 
close to 1 — a/2- Another higher commensurate phase is the 
(2 X 1) (Fig. |4(c)[ ) which is generated by a translation of the 
basis with the reciprocal lattice vectors of the a c(2 x 1) lattice 
H. Finally, the (2\/2 x \/2) phase (Fig. |4(e)l i is similar to 
the (2 X 1) phase. The lattice is generated by a translation 
of the basis with vectors which are rotated 45° with respect 
to the pinning potential and the magnitudes of the vectors are 
2\/2as and \/2as. The phase is favored for mismatch values 
close to -0.27. 

The transitions between the different phases are found by 
investigating the positions and the heights of the peaks in the 
structure factor The results of our extensive numerical cal- 
culations are summarized in the phase diagram of Fig. |5(a)| 
which has been taken from Ref . 13511 . 



1. Nonlinear Response of the (1 X 1) phase 



tive single particle behavior which is expected to describe the 
threshold behavior for a commensurate phase in absence of 
defects, is inde pen dent of the dimension, leading to a depin- 
ning exponent lil [H [Hill C = 1/2. We also find that 
the critical force increases with the pinning strength and for 
6m —0.3 has linear dependence on the pinning strength, 
while for 6^ S —0.3 its dependence on pinning becomes 
sub-linear. Finally, we also note that for large driving forces 
f ^ fc the system is totally depinned and the dependence of 
the drift velocity on the force follows Eq. (|2]i with a linear 
dependence. 

B. Phase Diagram for 2D System 

Next we present a summary of the static properties of the 
2D PFC model in the presence of the external pinning poten- 
tial y) = Vb(cos(fcsx) +cos(fcsy)). The parameters cho- 
sen are r = —1/4 and ij; — —1/4. Analytically, we consider 
the density to be a sum of hexagonal and square modes, 

ij{x,y) = At(cos(xV3/2)cos(y/2) - icos(2/)) (9) 

+ Asi{cos{ksx) + cos{ksy)) 
+ As2 cos(fcsx) cos{ksy) 

k k - 

+ Ac cos{-^x) cos{-^x) + tp. 

For small values of the pinning strength the system is in a 
hexagonal I phase for all mismatches (Fig. |4(a)| i. When the 
pinning strength is large enough, the system will be in one of 
the commensurate phases, of which the (1x1) phase is an 
exact match with the pinning potential here (Fig. |4(b)| i. 

The other ordered phases are higher commensurate phases, 
which exist only when one of the reciprocal lattice vectors 
for the commensurate phase is close to the wave vectors of 
the square pinning potential Issll . One of these phases is the 



We now turn to the influence of an external driving force for 
the pinned, commensurate (1 x 1) and c(2 x 2) phases. De- 
pending on the values of Sm and Vq, for both phases we find 
both discontinuous and continuous depinning transitions. For 
Sm > —0.2 and Vq > 0.09 both continuous and discontinu- 
ous depinning mechanisms were found for the commensurate 
(1x1) phase. For smaller values of the pinning strength close 
to the IC transition, only discontinuous depinning transitions 
were found. Similar to the ID case, we identify two values 
of the critical forces for a discontinuous depinning, namely 
for when the force is increased and /^"^ when the force 
is decreased back to zero. Figure [6(a)| shows the velocity de- 
pendence with respect to the applied force for a discontinuous 
transition. We have also tested the effect of thermal fluctua- 
tions here and find that for temperatures low enough the hys- 
teresis remains unchanged up to some value which depends on 
the mismatch and pinning strength. The gap A/c = — /^*^ 
for a given 6n-i decreases as the pinning strength is increased 
(Fig. |6(b)[ ). Finally, when the gap A/c vanishes the depinning 
transition becomes continuous. In this regime we find that the 
sliding velocity follows a power low Vd oc (/ — fc)'' (Fig. [T]), 
consistent with ( = 0.5 in all cases studied here as in the ID 
case. We note that for 6^ < —0.2 only continuous depinning 
transitions were found for aU values of the pinning strength. 

Both depinning mechanisms are accompanied by structural 
changes. The system changes from a commensurate (1x1) 
phase (below critical threshold) to a distorted hexagonal phase 
(Figs.[8l9]l. 

2. Nonlinear Response of the c(2 X 2) phase 

Similar to the (1 x 1) phase, the commensurate c(2 x 2) 
phase also exhibits both discontinuous and continuous depin- 
ning. For low value of pinning strength discontinuous depin- 
ning transition was found (Fig. [TOb . while for large values of 
the pinning strength the depinning becomes continuous (Fig. 
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FIG. 5; The phase diagram in terms of the pinning strength (Vo) and 
mismatch 5-^ calculated (a) numerically, according to Ref. (3^ and 
(b) analytically using approximation of the density given by Eq. l[9). 
The insets in (a) and (b) show the phase diagram close to 5^ = 0. 
The circles in (b) mark the values for which the approximation for 
the density given by Eq. ([9} breaks down. 




FIG. 4: The phases that minimize the free-energy, according to Ref. 
(3^ : (a) hexagonal, (b) square (1 x 1), (c) square (2 x 1), (d) square 
c(2 X 2), and (e) square 2V2 X \/2. The upper panels represent 
the density plotted in a gray colormap and the corresponding lattice 
vectors, while the lower panels show the structure factors and the 
relevant reciprocal lattice vectors. The black contours in Figs. |4(c)| 
and |4(e)| show the bases which generate the (2 x 1) and (2\/2 x \72) 
lattices. 



fTTT i with the exponent C, ~ 0.5. In both cases the structure of 
the systems changes when the force is applied from a C phase 
to a distorted hexagonal depinned phase (Figs. I12I13I ). 

While for the cases presented here the structural changes 
occur when the system starts to depin, it is also possible for 
changes to occur for forces below the critical threshold. For 
some values of the mismatch and pinning strength, if the force 
is rotated 45 degrees the c(2 x 2) phase will change first to a 
(1x1) phase from which the systems depins continuously 
as described before. Other force induced transitions between 
commensurate phases with no sliding are also possible. 

The hysteresis behavior of the depinning transition found 
for sufficiently small pinning strength Vb and the critical ex- 
ponent ^ for the continuous transition for larger V) are con- 
firmed by calculations of peak velocity Vp from Eq. (Q. For 
large Vq where there is no hysteresis, the behavior of Vp as 
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FIG. 6: (a) The variation of the velocity with respect to the external 
force for a discontinuous transition for the commensurate (1x1) 
phase (5m = 0.125, Vo — 0.0350) and (b) the variation of the gap 
A/c vs Vo for the same mismatch ~ 0.125. 



function of / show a depinnning transition at a critical force 
fc- A power-law fit of the velocity near fc gives an exponent 
( = 0.52 ± 0.03 which is consistent with the estimate using 
the velocity definition in Eq. Q. 

The determination of the velocity response from the peak 
positions allows us also to study the response to an additional 
force fy applied perpendicular to the longitudinal force in 
the moving state. For fx >> fc the longitudinal velocity Vx 
is proportional to the force since in the moving state the exter- 
nal pinning potential in the direction of the force appears as a 
time dependent perturbation in a reference system co-moving 
with the lattice, with a vanishing time average ll40ll . However, 
the pinning potential remains static in the transverse direc- 
tion ll40ll . One then expects that for sufficiently larger pinning 
strength, a transverse depinning transition is possible for in- 
creasing force fy while fx is kept fixed. Fig. \T4[a) shows the 
behavior of the transverse velocity component Vy when an in- 
creasing fy is applied in the moving state with fixed fx > fc- 
The transverse critical force fyc decreases with the longitudi- 
nal force fx and appears to vanishes at the longitudinal depin- 
ning transition fc, as shown in Fig. [TJb). 



FIG. 7: Dependence of the velocity on the external force for a con- 
tinuous transition for the (1 x 1) phase ((5,n = 0.125, Vo = 0.0900). 
The vertical arrow marks the critical force fc. The triangles represent 
the numerical data, while the continuous line is a power-law fit with 
C = 0.50 ±0.03. 



IV. DISCUSSION AND CONCLUSIONS 

In this work we have considered the recently developed 
phase field crystal model |f32| l in the presence of an external 
periodic pinning potential [34,3^ and a driving force. As the 
model naturally incorporates both elastic and plastic deforma- 
tions, it provides a continuum description of lattice systems 
such as adsorbed atomic layers on surfaces or 2D vortex lat- 
tices in superconducting thin films, while still retaining the 
discrete lattice symmetry of the solid phase. The main advan- 
tage of the model as compared to traditional approaches is that 
despite retaining spatial resolution at its lowest length scale 
its temporal evolution naturally follows diffusive time scales. 
Thus the numerical simulation studies of the dynamics of the 
systems can be achieved over realistic time scales, which for 
example in the case of adsorbed atomic systems may coitc- 
spond to many orders of magnitude over the time scale used 
in microscopic atomic models. In this work we have exploited 
this method to determine the phase diagram in one and two di- 
mensions as a function of lattice mismatch, pinning strength 
and a driving force. 

We have concentrated on the nonlinear response to an ex- 
ternal driving force on the most common stable commensu- 
rate states, namely the (1 x 1) and the c(2 x 2) phases. These 
are particularly interesting cases which are relevant for phys- 
ical systems of current interest and accessible experimentally 
such as, driven adsorbed layers ||^|2^, which determine the 
sliding friction behavior between two surfaces with a lubricant 
I IC I and between adsorbed layers and an oscillating substrate 

Our results for the phase field crystal model with over- 
damped dynamics indicate both discontinuous and continu- 
ous transitions depending on the magnitude of the pinning 
strength. For high enough pinning strengths continuous tran- 
sitions occurred with the velocity near the transition scaling as 
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FIG. 8: Change of the lattice structure (upper panels) and the cor- 
responding structure factor (lower panels) with the applied force for 
5m = 0.125, Vb ~ 0.0350, where depinning is discontinuous. Im- 
age (b) corresponds to / = 0.11 with a non-moving initial config- 
uration, while for (c) the applied force is the same but the initial 
configuration is a moving one. The case (a) / — 0.07, (d) / = 0.13 
are outside of the hysteresis region and same result is obtained with 
moving or non-moving initial configuration. 



(/ — /c)^/^, independent of the dimension of the system. This 
is as expected, since for a commensurate state in a strong peri- 
odic pinning potential, each 'particle' acts independently and 
the model reduces to an effective single particle in a periodic 
potential, with a known depinning exponent of 1/2. Perhaps 
more interesting is the observation of discrete transitions and 
hysteresis loops found at low pinning strengths. In the two- 
dimensional case, the observed hysteresis behavior is consis- 
tent with the arguments and atomistic molecular dynamics 
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FIG. 9: Change of the structure factor with the applied force for 
5m ~ 0.125, Vb — 0.0900, where depinning is continuous. The im- 
ages correspond to (a) / = 0.3776 (before the depinning transition 
marked by the vertical arrow in Fig. |7(a)[ (. (b) / — 0.3782 (after the 
depinning transition). 



simulations of driven adsorbed layers l26ll indicating that 
hysteresis remains in the overdamped limit. However, our re- 
sults show that it disappears for large enough pinning strength. 
For the discontinuous transition, there are two different criti- 
cal values, f"'^ > Z^*^, correspond to the static and kinetic crit- 
ical forces, respectively, which lead the to stick and slip mo- 
tion at low sliding velocities as observed experimentally ifioll . 
The general features observed of hysteresis and power-law be- 
havior near the continuous depinning transition are also of in- 
terest for driven charge density waves lITsl l29ll and driven flux 
lattices ll4lll43ll although in these cases there are important ad- 
ditional effects due to a large degree of disorder in the pinning 
potential. Whether the present phase field model also shows 
the same main features observed for the atomistic model in 
presence of thermal fluctuations ll26l l27l l28ll is an interesting 
question which will require further investigation. 

In addition to the longitudinal depinning transition where 
the lattice system is moving in the same direction as the driv- 
ing force, a driven two-dimensional lattice on a periodic po- 
tential can also show an interesting behavior for the trans- 
verse response in the moving state. When the lattice is al- 
ready moving along some symmetry direction of the pinning 
potential the response to an additional force applied in the 
direction perpendicular to the longitudinal driving force may 
lead to a depinning transition for increasing transverse force 
ll40ll . Such transverse depinning has been found in different 
driven lattice systems with periodic pinning including driven 
vortex lattices iiH l42ll and adsorbed layers |27J in standard 
molecular dynamics simulations. In the present PFC model, 
we have obtained similar results for the transverse depin- 
ning. Experimentally, some evidence of transverse pinning 
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FIG. 10: (a) The variation of the velocity with respect to the external 
force for a discontinuou.s transition for the c(2 x 2) phase ((5m — 
-0.50, Vb = 0.099) and (b) Afc vs Vo for 5™ = -0.50. 

has been observed in measurements on charge-density waves 
iH . Wigner solid iH] and vortex lattices 14^, although in 
these cases disorder in the pinning potential plays a more im- 
portant role. 

While the hysteresis behavior is expected in the presence 
of inertial terms both in ID and 2D, it is quite interesting to 
see it in ID when the dynamics being used is overdamped and 
purely relaxational. In 2D, it can be argued that topological 
defects such as dislocations can lead to this behavior even with 



overdamped dynamics but these defects are not available in 
ID. It is interesting to speculate that the hysteresis behavior 
is intimately related to the need for plastic deformations to 
mediate the transition from one lattice structure to another 
Work on these problems is already in progress. 
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